Comment on "On two-dimensional magnetohydrodynamic turbulence" [Phys. 

Plasmas, 8, 3282 (2001)] 
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Biskamp and Schwarz [Phys. Plasmas, 8, 3282 (2001)] have reported that the energy spectrum 
of two-dimensional magnetohydrodynamic turbulence is proportional to k~^^^, which is a predic- 
tion of Iroshnikov-Kraichnan phenomenology. In this comment we report some earlier results which 
conclusively show that for two-dimensional magnetohydrodynamic turbulence, Kolmogorov-like phe- 
nomenology (spectral index 5/3) is better model than Iroshnikov-Kraichnan phenomenology; these 
results are based on energy flux analysis. 



In a recent paper, Biskamp and Schwarz (now on 
referred to as Ref. 1) discuss energy spectrum and struc- 
ture functions for two-dimensional (2D) magnetohydro- 
dynamic (MHD) turbulence. Based on numerical calcula- 
tion of energy spectrum, they claim that the spectrum of 
2D MHD turbulence agrees with Iroshnikov-Kraichnan's 
(IK) fc^'^/^ law with some modifications. The purpose 
of the present comment is to show that the above claim 
is inconclusive. We bring to notice here an alternative 
point of view based on energy cascade rates which sup- 
ports Kolmogorov's energy spectrum K~^^^ for 2D as 
well as 3D MHD turbulence; this result is reported in 
Verma et al. § and Dar ||. 

MHD turbulence phenomenologies are discussed in 
Verma et al. (referred to as Ref. 2) . For zero velocity- 
magnetic correlation (u ■ b) , IK phenomenology predicts 



i;"(fc) = E^k) = A(HVSi)^/2fc"^/^ 
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where H is the total energy flux, Va is the Alfven veloc- 
ity, and A is a universal constant. Dobrowolny et al. 
have generalized IK's arguments for nonzero u • b and 
showed that the energy cascade rates H^ of = u ± b 
are equal irrespective of and E~ ratio, i.e.. 
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Marsch g proposed Kolmogorov-like phenomenology in 
which 
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This is a also a limiting case of a generalized phenomenol- 
ogy of Matthaeus and Zhou H . Clearly, 
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Verma et al. [||, Dar and 
have numerically computed the 
spectral exponents for 2D MHD turbulence. Biskamp 
and Welter |Ql support IK's k'^^'^ energy spectra, but 
Verma et al. ^ and Dar find numerical uncertainties 
too significant to be able to distinguish between the ex- 
ponents 3/2 and 5/3 (see Fig. 1 of Ref . 2). Biskamp and 
Schwarz |ll| do not provide the error bars for the spectral 
indices (see Fig. 7 of Ref. 1). Since 5/3 and 3/2 are so 
close, the claims of Biskamp and Schwarz in favor of 3/2 
may not be conclusive. As stated by them, the intermit- 
tency exponents do not clarify the matter any further. 
On the other hand, based on energy flux studies, Verma 
et al. and Dar could show quite conclusively that 
Kolmogorov-like phenomenology models 2D MHD turbu- 
lence better than IK phenomenology. 

Verma et al. ||] and Dar |^ numerically computed 
the energy fluxes n+ and H" for various E~/E^ ra- 
tios. The cascade rates of majority species (larger of E^ 
and E~^) was always found to be greater that those of 
minority species. To illustrate we have plotted H^ for 
E-/E+ w 0.2 in Fig. 1 (taken from Dar §). The same 
results are observed in Verma et al. however, the er- 
ror bars in Dar |^ is relatively smaller (w 5%) because of 
better averaging. Clearly n+ > H^. Incongruence of the 
above result with the IK predictions [Eq. (2)] clearly in- 
dicated that IK phenomenology is not valid for 2D MHD 
turbulence. For E~ / E'^ in the range of 0.2 to 1, Verma 
et al. 01 and Dar 01 find that 



E- 
E+ 



(5) 



This result is in agreement with the predictions of 
Kolmogorov-like phenomenology for MHD turbulence 
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with = K~ [Eq. (4)]. Intermittency corrections to 
energy fluxes are typically small. Hence, we do not ex- 
pect the large difference in energy fluxes to result from 
intermittency. This way Verma et al. and Dar 
showed that for 2D MHD turbulence, Kolmogorov-like 
phenomenology is better model than IK phenomenology. 
Recent advances on theoretical fronts also tend to indi- 
cate that Eq. (4) may even be valid for smaller E~ /E^ 
with K+ ^K-. 

Current theoretical and numerical papers argue 
that Kolmogorov's energy spectrum in MHD is due to 
local Alfven effects. The Alfven waves are scattered by 
the "local mean magnetic field", rather than the global 
mean magnetic field. Hence, the effective time-scale will 
be comparable to the nonlocal time scale resulting in Kol- 



mogorov's energy spectrum for MHD turbulence. The 
above argument is expected to hold in both 2D and 3D. 

To conclude, Verma et al. ||] and Dar's |^ results based 
on energy fluxes support Kolmogorov's spectrum for 2D 
MHD turbulence. We believe Biskamp and Schwarz's 
claim favoring k"^^"^ energy spectrum is incorrect. 

Note added after the receipt of the Reply: The authors 
thank Prof. Biskamp for pointing out Grappin et a/.'s 
paper |10|. However, for high normalized cross helicity 
(w 0.9), Verma et al. and Dar find the exponents 
of « fc^™ in the range of 1.5-1.7, but the ratio 
H+/H^ ;2> 1 (5 to 10) . Hence numerical results of Verma 
et al. Pi an d Dar |^ are not in agreement with Grappin 
et al. 'spc|] predictions that H+/H^ w m'^ /m^ . 



D. Biskamp and E. Schwarz, Phys. Plasma 8, 3282 (2001). 
M. K. Verma et al, J. Geophys. Res. 101, 21619 (1996). 
G. Dar, Ph.D. thesis, I. I. T. Kanpur, 2000. 

M. Dobrowlny, A. Mangeney, and P. Veltri, Phys. Rev. Lett. 45, 144 (1980). 

E. Marsch, in Reviews in Modem Astronomy, edited by G. Klare (Springer- Verlag, Berlin, 1990), p. 43. 
W. H. Matthaeus and Y. Zhou, Phys. Fluids B 1, 1929 (1989). 

D. Biskamp and H. Welter, Phys. Fluids B 1, 1964 (1989). 
M. K. Verma, Phys. Plasma 6, 1455 (1999)^ 



J. Cho, A. Lazarian, and E. T. Vishniac, astro-ph/0105235 (accepted in Astrophys. J.) 



R. Grappin, A. Pouquet, and J. Leorat, Astron. Astrophys. 126, 51 (1983). 




-0.01 ' ' ' — ' — ' — ' ' — ' — ' — 

1 10 100 

k 



FIG. 1. Energy fluxes IT versus k for a two-dimensional run with Bo = 0.0 and E~/E+ = 0.18. Here n+ > H". The figure 
is taken from Dar [3]. 
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